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For  several  years  RAND  has  Investigated  meens  of  applying  principles  of 
mass-transfer  and  ablation  coollzig  to  problems  of  atmospheric  re-entry 
and  to  the  design  of  efficient  hypersonic  flight  vehicles.  In  connection 
with  ablation  cooling,  this  research  memorandum  examines  the  binary 
boundary-layer  problem.  The  study  shoiild  make  possible  more  rapid  and 
reliable  estimates  of  surface -cooling  methods  for  use  with  hypersonic 
vehicles,  euoh  as  Intercontinental  balllstlo  missile  nose  cones. 

A  binary  boundary  layer  Is  one  in  which  some  foreign  substance  has  been 
Injected  to  alter  the  properties  of  the  flow,  notably  Its  heat-transfer 
oharacterlstlcs.  Several  methods  available  for  acconpllshlng  this 
Injection,  or  mass  transfer,  are  the  transpiration  of  gas  through  slots, 
emd  the  ablation  or  sloughing  away  of  surface  particles. 

The  meohanl.em  of  laminar  binary  boundary  layer  flow  is  discussed  In 
mathematical  terms,  and  five  different  analyses  Involving  a  variety  of 
injected  substances  are  reviewed.  Generalized  e^giresslons  are  then  developed 
for  predicting  heat-transfer  and  skln-frlctlon  performance  In  the  presence 
of  mass-transfer  cooling  for  laminar  flow  over  a  flat  plate.  The  results 
indicate  that  different  foreign  materials  (for  example,  hydrogen,  carbon 
dioxide,  and  iodine  vapor)  Injected  into  the  boundazy-layer  stream  reduce 
heat-transfer  and  skln-frlctlon  coefficients  by  an  amount  idilch  depends  on 
the  molecular  weight  of  the  Injected  material.  In  conclusion,  mass-transfer 
cooling  In  a  turbulent  boundary  layer  and  sublimation  cooling  are  considered. 
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SUMMARY 

The  purpose  of  this  report  Is  to  examine  the  several  theoretical 
analyses  of  laminar  and  ttirhulent  binary  boundary  layers.  A  generalization 
of  the  different  theories  leads  to  sio^lifled  expressions  for  the  heat 
transfer  and  friction  In  both  laminar  and  turbulent  binary  boundary 
layers.  The  results  Indicate  that  the  effects  of  different  foreign 
materials  are  primarily  dependent  on  their  relative  molecular  weights . 

The  importance  of  the  Injection  of  foreign  materials  to  the  stability 
eharaoterlstlcs  of  the  flow  is  discussed  and  experimental  results  seem  to 
show  non-transltlonal  flow  for  moderate  Injection  rates  even  with  light 
gases . 

Sublimation  cooling  Is  treated  and  typical  numerical  results  are 
presented  for  three  materials «  ice,  carbon  dioxide,  and  ferrous  chloride. 
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SYMBOLS 

Chapman-Rubealn  parameter  defined  In  Eq,  (17) 
local  skin  friction  coefficient 

local  skin  friction  evaluated  for  solid  vnsill  exposed  to  some 
free  stream  conditions  and  held  at  same  temperature  as  the 
actvial  wall 

local  heat-transfer  Stanton  number,  q/p^u  C  (T  -  T  ) 

^eep  r  w 
e 

local  Stanton  number  evaluated  for  solid  wall  exposed  to  same 
free  stream  conditions  and  held  at  same  teniperature  as  the 
actuial  wall 

mass -transfer  Stanton  number,  in/ftn  Y. 
specific  heat  at  constant  pressure 
specific  heat  at  constant  volume 
ordinary  diffusion  coefficient 
dimensionless  stream  function 
gravitational  constant,  52.2 
enthalpy  associated  with  change  of  phase 
heat-transfer  coefficient.  Eg.  (19) 
thermal  conductivity 

refers  to  wedge  flows,  where  free  stream  velocity,  u  ,  varies 
as  x“  * 

molecule  weight  of  coolant  gas 

molecule  weight  of  pure  air 

ratio  of  molecule  weights,  m^/ffig 

mass  flow  rate  of  coolant  gas  at  position  x  along  the  surface 
pressuire 

partieil  pressure  of  coolant  gas 
local  heat-transfer  rate  per  unit  area, 
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q  loco,l  heat -transfer  rate  per  \mlt  area  evaluated  for  solid  vail 

exposed  to  seune  free  stream  conditions  and  held  at  same 
temperature  as  actual  wall 

R  universal  gas  constant 

r  recovery  factor 

T  teu^erature 

recovery  temperature,  defined  as  the  wall  tanmerature  where 

iV(aT/ay)^  -  0 

T’l'  reference  temperature  defined  in  Eq,  (9) 
u  eon^'  nent  of  velocity  parallel  to  surface 

defined  In  Eq.  (23) 

V  component  of  velocity  normal  to  surface 

X  coordinates  along  the  body 

V  mass  fraction  of  coolant  gas 

y  coordinates  normEd  to  the  body 

Greek  Symbols 

a  condensation  coefficient 

y  ratio  of  specific  heats 

b*  displacement  thickness,  Eq.  (22) 

irj  transformed  coordinate  defined  in  text 

6  dimensionless  temperature 

\  defined  in  Eq,  (17) 

dynamic  viscosity 

V  kinematic  viscosity 

p  density 

^  dimension  mass  fraction 


stream  fvinotion 
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PlmenslonleBa  Nuaibers 

M  Mach  number;  ratio  of  local  speed  to  local  speed  of  sound 

Pr  Prandtl  number,  uc  /k 

r  "I  ^ 

^®inln  critical  Reynolds  number  based  on  u^ 

Re„  local  Reynolds  number,  u  x/v^ 

Sc  Schmidt  ntimber,  v  /D^ 

Subscripts 

1  refers  to  pure  coolant 

2  refers  to  pure  air 

e  evaluated  at  outer  edge  of  the  boundary  layer 
r  refers  to  recovery  conditions,  i.e.,  vhere  k(dT/dy)^  ■  0 
V  evaluated  at  vail  conditions 

X  refers  to  local  conditions 

y^  evaluated  at  y^  as  defined  in  Eqt  (23) 

Superscript 


*  evalmted  at  the  reference  temperature,  Eq.  (9} 
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I.  INTRODUCTION 


Recently,  the  alleviation  of  the  high  heating  rates  encountered  hy 
surfaces  of  hypersonic  vehicles  has  been  recognized  as  an  Important 
problem.  One  of  the  cooling  methods  that  appears  to  have  great  ultimate 
promise  is  mass -transfer  cooling,  wherein  a  "foreign"  material  Is 
transferred  from  the  surface  Into  the  boundary  layer.  This  has  a  two¬ 
fold  advantage  in  alleviation  of  the  heat-transfer  problem.  The  trans¬ 
ferred  coolant  may  absorb  heat  from  the  boundary  layer  throu£^  a  phase 
oheuQge  (sublimation,  evaporation,  melting,  etc.)  and/or  by  acting  as  a 
dispersed  heat  sink.  Zt  will  be  advantageous,  therefore,  to  employ 
coolants  with  hl£^  heats  of  svCbllmatlon  (or  evaporation,  melting,  etc.) 
as  veil  as  high  thermal  heat  capacities.  In  addition,  it  has  been 
shown  that  the  Introduction  of  a  material  (with  Its  normal  velocity 
oong>onent)  at  the  surface  acts  to  decelerate  the  flew  and  consequently, 
to  reduce  the  skin  friction.  This  also  Inplles  a  reduction  In  heat 
transfer  at  the  wall. 

The  usual  boundary  layer  eq\iatlons  are  complicated  by  the  appearance 
of  (a)  an  equation  defining  the  conservation  of  the  species  at  any  point 
in  the  boundary  layer,  and  (b)  transport  terms  which  result  from  thenso- 
dynamlo  coipllng  coefficients  such  as  the  thermal  diffusion  coefficient. 

There  are  a  number  of  methods  for  effecting  the  Inje-tlon  of  a 
foreign  material  Into  the  boundary  layer,  and  the  following  descriptions 
have  been  advanced  to  describe  specific  mass-transfer  cooling  schemes s 
1.  Transpiration  cooling 
8.  Film  cooling 

a)  using  liquid  as  a  coolant 
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b)  uBlng  gas  as  a  coolant 

3 .  Ablation 

a )  sublimation 

b)  other  ablation  phenomena,  such  as  melting,  erosion,  ftiflion, 
etc. 

(l)  * 

These  schemes  are  diagramatleally  shown  in  Fig.  1.  With  the 

exception  of  film  cooling  with  a  gas,  these  methods  all  Involve  the  same 
mechanism  in  the  gaseous  phase  of  the  boundary  layer.  However,  there  is 
a  difference  In  the  boundary  conditions  at  the  wall  which  distinguishes 
the  three  methods  from  a  thermodynamic  as  well  as  a  mechanical  viewpoint. 
Transpiration  cooling  Involves  the  introduction  of  a  coolant  gas  through 
a  pororiB  surface.  Consequently,  the  rate  of  fluid  injection  through  the 
surface  and  into  the  boundeory  layer  may  be  arbitrarily  adjusted  by 
purely  mechanical  means  and  the  temperature  at  the  surface  may  thereby 
be  regulated  depending  upon  the  injection  rate  and  storage  temperature  of 
the  coolant  gas.  It  should  be  noted  that  a  transpiration-cooling  system 
requires  pvmps,  storage  tanks,  pressure  regulators,  and  accessory 
plumbing.  In  addition,  the  fabrication  and  maintenance  of  porous  surfaces 
represent  a  difficult  engineering  problem. 

A  film-cooling  system  involves  punping  a  liquid  or  gas  onto  the 
surface  through  an  inlet  slot  configuration  such  that  a  thin  film  of 
the  materled  covers  the  surface.  This  acts  as  an  insulating  coating  and, 
in  the  case  of  a  liquid,  absorbs  heat  in  its  vaporization.  These  systems 
are  usually  limited  by  such  characteristics  as  the  stability  of  the 

figures  in  parentheses  refer  to  references  presented  at  the  end  of 
the  text. 
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Various  cooling  methods 
Hot  gas - ►  Hot  gas - ► 

^  ^  ^  ^ 


’4  !//////  4 

Coolant 
Film  cooling 
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Coolant 

T ronspirotion  cooling 


Liquid  film  cooling 


Hot  gas 


Hot  gos 
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Solid  Solid 

Ablation  cooling 

Fig.  I  —  Various  cooling  methods 
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liquid  film  and  the  pumping  power  available.  Fllm-coollng  systems 
reqtdre  essentially  all  the  plmblng  and  control  equipment  of  a  trans- 
plratlon-cooled  operation}  however^  the  surface  construction  Is 
mechanically  simpler. 

A  BUbllmatlon-aoollng  system  Is  self -controlled  through  the 
relation  between  the  vapor  pressure  and  surface  teaperature  of  a  solid, 
bhe  Clauslus-Clapeyron  equation! 

In  +  B  (l) 

where 

p^  •  vapor  pressure  of  subliming  material 
R  ■  universal  gas  constant 
R.  ■  heat  of  atR)llfflatl6n  per  mole 
B  ■  constant  of  Integration 

Beat  Is  absorbed  by  the  material  as  It  sublimes.  Thus,  the  heat  transfer 
Into  the  Interior  Is  reduced  in  two  ways:  (a)  direct  absorption  In  tbs 
form  of  heat  of  sublimation,  and  (b)  reduction  of  heat  transfer  because 
of  the  movement  of  the  sublimed  mass  away  from  the  surface.  The  mass 
release  at  the  svirface  depend^  \pon  the  heat  of  sublimation  and  the 
temperature  of  the  stxrface.  Furthermore,  the  surface  temperature  can  no 
longer  be  arbitrarily  controlled  and.  In  fact,  will  always  find  ''its  own 
level"  depending  upon  the  heat  load,  heat  of  sublimation,  and  the 
ejcternal  flow  situation.  It  should  be  noted  that  film  cooling  with  a 
liquid  Is  essentially  a  sublimation  process  provided  that  the  surface  Is 
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completely  covered  by  a  liquid  layer. 

In  addition  to  sublimation,  more  complex  ablation  cooling  Bchemee 
may  be  viflualizecL.  Depending  upon  the  surface  material  emd  the  fli^t 
cotidltions,  it  is  possible  to  have  such  phenomena  as  fusion  of  the 
sxtrface  material,  mechanical  erosion,  diaseclatlen  of  both  air  and 
surface  material,  ionization,  and  chemical  reactions  between  the  oospoo- 
nents  in  the  boundary  layer.  The  obvious  conplicatlons  Involved  in  an 
analysis  of  these  cooplex  ablating  systems  have  prevented  any  really 
accurate  description  of  the  mechanism. 

There  is  another  method  of  classifying  these  systems  which  may  be 
helpful  conceptually;  i.e.,  by  specifying  the  method  of  controlling  the 
rate  of  injection.  In  the  ease  of  transpiration  cooling,  as  we  have 
seen,  the  rate  is  arbitrarily  controlled  depending  upon  certain  mechanical 
requirements  such  as  the  porous  surface  and  pumping  power  available.  This 
may  be  considered  an  arbitrarily  controlled  system.  For  a  subliming  or 
ablating  surface,  however,  the  rate  of  injection  Is  determined  by  the 
heat  of  sublimation  and  the  siirface  temperature.  For  a  set  of  flight 
oondltlons  and  a  surface  material,  the  steady-state  Injection  Is  fixed 
thermodynamically.  We  may  call  this  a  self -controlled  system.  Film 
cooling  may  he  considered  a  mixed  system,  with  the  rate  of  mass  transfer 
into  the  boundary  layer  being  self -controlled,  but  the  flow  rate  of 
liquid  or  gas  over  the  surface  remaining  arbitrary. 
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II.  BINARY  LAMINAR  BOUNDARY  I-AYER  THEORY 

Baron/'^'^^  Eckert,  et  al. Szlklaa  and  Banas,^^^  and  Qroas^®^ 
have  Investigated  theoretically  the  problem  of  mass -transfer  cooling  in 
the  laminar  boundary  layer  on  a  flat  plate.  The  introduction  of  a 
species  conservation  equation  as  well  as  the  appearemce  of  thermodynamic 
covgjllng  terms  in  both  the  species  and  energy  equations  complicates  the 
mathematical  analysis  of  the  boundary  layer. 

(o) 

The  equations  have  been  derived  by  Hall'^'  who  first  treats  the 
multicon5)onent  fluid  system  and  then  from  an  order -of -magnitude  argument 
obtains  the  boundary  layer  equations  for  a  binary  boundary  layer.  The 
final  binary  boundary  layer  equations  for  flow  over  a  flat  plate,  neg¬ 
lecting  the  effects  of  thermal  diffusion  may  be  reproduced  as  follows: 


continuity:  (pu)  (pv)  «  0  (a) 

Momentum:  fy  “  ^  3y^ 

Energy: 

f  +  P  D^2(Cpi  -  Cpg)  |I  (4) 

Species:  pu  0  +  pv  g  «  ip  Tp  |l  1  (5 ) 


The  botuidary  conditions  for  this  system  of  equations  follow  from 
physical  considerations: 
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u  =  o 

T  =  T 
w 

Y  =  Y 

w 

u  «  u 

e 

T«T>  aty— 00  X^) 

Y  «a  0 

Since  we  are  concerned  with  a  aeventh>order  syetem  of  equations^ 
another  boundary  condition  mtist  be  specified.  This  condition  nay  be 
obtained  by  noting  that  the  nass  velocity  of  the  air  molecules  disappears 
at  the  surface  of  the  plate)  that  Is  to  say,  there  is  no  net  mass  transfer 
of  the  boxmdary  layer  air  Into  the  plate  surface.  Therefore,  the  nass 
flow  of  air  by  convection  away  from  the  surface  must  be  equal  and 
opposite  to  the  diffusive  flow  of  air  toward  the  surface.  This  consid¬ 
eration  yields  the  following  bovindary  condition: 

V  =  -  y  ^  at  y  o  0  (8) 

This  system  of  Eqs.  (2)  -  (8),  forms  the  starting  point  for  the 
various  investigators  cited  above.  In  all  cases,  a  tranoforraation  of 
the  coordinates  is  next  introduced  with  the  result  that  the  system  of 
partial  differential  equations  Is  changed  Into  a  new  set  of  Interdependent 
ordinary  differential  equations.  This  system  of  equations  is  still 
difficxdt  to  solve  since,  in  general,  all  of  the  physical  properties  are 
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functions  of  the  local  temperat\ira  and  concentration  of  the  particular 
gas  mixtures  being  investigated.  To  obtain  a  representative  nxnnber  of 
solutions  in  a  reasonable  time  requires  the  use  of  high-speed  electronic 
oonputers.  Although  the  some  basic  system  of  differential  equations  was 
used  by  the  various  investigators (Eqs.  (2)  -  (8)),  somewhat  different 
aesuiqptlons  wei's  imposed  to  obtain  the  final  solution.  A  brief  review 
of  the  available  analyses  is  given  below. 
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III.  RSVIEW  OF  AVAHABIaE  BINARY  LAMINAR  BOUHDARY  LAYER  ANALYSES 
CON8TAMT  PROPERmf  MALYSIS^^^^ 

Considerable  Insight  Into  the  mass -transfer  cooling  process  Is 
obtained  If  It  Is  assimed  that  the  Injected  coolant  has  physical 
properties  not  markedly  different  from  those  of  the  main  stream,  thereby 
permitting  the  assvoqptlon  of  constant  physical  properties.  Ilie  advantage 
of  such  a  constant  property  solution  Is  that  the  dimensionless  heat- 
transfer  and  ekln-frlotlon  results  ere  dependent  on  a  mlnlnun  number  of 
parameters.  This  Is  Indicated  in  Figs.  2  and  3  vhloh  ooopara  the 
dimensionless  quantities  Vhlch  are  of  Importance  for  the  solld-vall, 
heat-transfer  case  and  those  arising  in  the  presence  of  mass  transfer. 

For  the  solid  vail.  It  has  been  demonstrated^^^  that  the  oonstant 
property  solutions  for  dimensionless  skin  friction,  Nusselt  numbers,  and 
recovery  factors  have  additional  value  In  that  they  may  be  used  even 
when  large  variations  In  physical  properties  are  encountered  (including 
dissociation),  provided  the  properties  are  evaluated  at  a  so-csdled 
reference  tenperature,  T*,  which  may  be  given  explicitly  In  terms  of  the 
surface,  the  free  stream,  and  the  recovery  temperatures : 

=  T  +  0.5  (T  -  T  )  +  0.22  (T  -  T  )  (q) 

e  '  V  e  re  ^ 

It  will  be  demonstrated  that  this  temperature,  1*,  will  also  prove  of 
value  In  correlating  the  mass -transfer  results. 

It  should  be  noted  from  Fig.  3  that  the  heat  transfer,  q,  is  defined 
In  terms  of  the  temperature  gradient  at  the  wall.  This  definition  Is 
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FOR  GIVEN  FREE  STREAM  CONDITIONS 
Tw=Tr  WHEN-k„(iL)^  =  0 

Tr“Tr(Me,Te) 


DEFINITION  OF  HEAT-TRANSFER  COEF FICIENT ,  h 

q«h  (Tw~Tr) 

Nu  s  — 7“ 
k 

-p=  -  f  (Me,Te,Tw) 

vRe^ 


CONSTANT  PROPERTIES, 


Nu% 


f  (Pr) 


Tr  =  Tr(Pi-) 

Fig.  2 — Flat  plate,  solid  wall,  no  moss  transfer 
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q 


FOR  GIVEN  FREE  STREAM  CONDITIONS 
T«.Tr  .0 

Tf'Tf  (Meje,  INJECTED  6As) 

DEFINITION  OF  HEAT-TRANSFER  COE FFICIENT ,  h 

q  ®  h  (Tw  -  Tf) 
hx 


Nu  * 


Nux 

-s/Rfix 


*  f  (l^e  ,Te,Ty^,  INJECTED  GAs) 

^  6  0 

Nu 


CONSTANT  PROPERTIES,  -4=  *  ^  ,  Pr  ) 


Tr=Tr(^ 


Fig.  3 — Flat  plate,  mass  transfer 
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oonvenient  in  that  this  q['uauitlty>  q,  represents  the  convective  heat 
transferred  to  the  surface  from  the  boundary  layer >  and  In  this  sense  the 
boundsry  layer  heat  transfer  is  considered  separately  from  the  enthalpy 
carried  across  the  surface  by  the  coolant.  To  be  consistent  with  this 
point  of  view,  the  recovery  teinperature,  Is  defined  as  that  tenqperature 
where  the  wall -temperature  gradient  vanishes}  In  this  case  there  still 
exists  a  transport  of  enthalpy  across  the  surface,  by  virtue  of  the 
coolant  flow,  but  there  Is  no  convective  heat  transferred  to  the  wall 
from  the  boiuidary  layer. 

Returning  to  the  solutions  of  Eqs.  (S)  -  (8)  It  is  obvious  that 
the  energy  equation  Is  linear  and  consequently,  the  general  solution  of 
the  cooplete  equation  may  be  obtained  from  the  addition  of  (a)  a  general 
solttblon  of  the  homogeneous  equation  (that  Is,  neglecting  the  dissipation 
term)  and  (b)  a  particular  solution  of  the  complete  equation.  The  par¬ 
ticular  solution  results  In  the  specification  of  the  recovery  factor,  a 
direct  measure  of  the  temperature  assumed  by  the  surface  If  It  Is 
allowed  to  come  to  equilibrium  with  the  surroundings  by  convection  alone. 

We  need,  therefore,  only  to  dlreot  our  attention  to  the  general  solution 
of  the  homogeneous  eq'uatlon.  To  accooplleh  this  solution,  a  stream 
function,  tjr  ,  Is  first  Introduced  to  satisfy  the  continuity  equation,  and 
a  new  Independent  variable,  77  ,  and  the  new  dependent  variables  defined 
below  are  substituted  Into  the  original  eqviatlons: 


u  = 


ay  ’ 


V  e 


ax 


stream  function: 


(10) 
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f  - 

7f  m  (y/x)yu^  x/v 
0  -  (T  -  \)/(T^  -  Ty) 
»$  -  (Y^  -  Y)Ay 


(11) 


The  following  equations  result* 


Momentum* 

■ 

1 

O 

(12) 

Energy* 

afo  f  0 

(13) 

Diffusion* 

a  <(  ^  1  qa  f  ^  ■  0 

(14) 

Boundary  conditions* 
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3  «  1 

^  ■  1 

An  iniportant  observation  eoosnon  to  all  flat -plate  binary  laminar 
bovindary  layer  solutions  is  that  the  mass  transfer  Into  the  boundary 
layer  must  vary  as  l/yiT  If  we  are  to  arrive  at  a  system  of  ordinary 
differential  equations.  Further^  we  have  assumed  an  Isothermal  surfaoe 
and  It  will  be  shown  that  this  la  conpletely  cootpatible  with  the  Imposed 
mass-transfer  distribution. 

The  velocity  profiles  for  the  constant  property  mass-transfer  system 
are  shown  In  Fig.  4  for  several  different  Injection  rates.  Inspection 
of  these  profiles  brings  out  the  following  conclusions t 

1.  The  effect  of  mass  addition  Is  to  thicken  the  velocity  boundary 
layer. 

2.  The  velocity  profile  beccmes  S-shaped  with  mass  addltlon>  and 
since  this  Is  known  to  be  an  vmstable  type  of  profile;  it  may  be  concluded 
that  mass  transfer  is  destabilizing. 

5.  The  boundary  layer  "lifts  off"  the  wall  at  a  relatively  low 
value  of  mass  transfer j  i.e.,  at  ^  ’^^u^  x/v  »  0.619 

Apparently  the  boundary  layer  equations  fall  to  describe  the  flow  field 
at  this  mass -transfer  condition* 

The  skin  friction  coefficient  and  Niisselt  nuniber;  presented  In 
Figs.  3  and  6;  are  seen  to  decrease  with  Increasing  mass  transfer;  both 
going  to  zero  at  the  limiting  value  where  the  boundary  layer  "leaves" 
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Fig. 4  —  Effect  of  air  injection  on  velocity  profile 
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Fig.  6 — Effect  of  air  injection  on  iocol  Nusselt  number 
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the  vail,  fflie  recovery  factor,  ahovn  in  Flg«  7>  la  aomewhat  reduced  by 
nasB  transfer  but  not  aa  markedly  as  the  Nusselt  values.  At  hyper sonio 
velocities  the  actual  heat  transfer  to  the  surface  Is  proportional  to  the 
product  of  the  Nusselt  nuD^ber  and  the  recovery  temperature,  and  ve 
conclude  that  a  considerable  reduction  In  heat  transfer  Is  obtainable 
with  modest  aoounts  of  coolant.  This  is  the  feature  vhleb  has  dravn 
attention  to  this  cooling  scheme. 

Xt  nay  be  noted  that  the  heat-transfer  solutions  reveal  that  the 
local  heat  rate,  is  proportional  to  l/y/%  ,  vhlch  is  precisely  the 
distribution  of  the  Injected  coolant.  A  slsple  heat  balance  on  a 
sublimation  or  transpiration  system  yields  the  following  expressloni 

\  K  •  "i)]  -  Sc 


where 


m  mass  flow  rate  of  coolant 
AH  ■  change  In  enthsilpy  due  to  phase  change 
c  «  specific  heat  of  the  pure  coolant 


factor 


RM-S516 

82 


Thus  for  the  eituation  where  the  coolant  enters  the  wall  at  a  constant 
temperature,  T^,  or  If  the  storface  la  subltalng,  the  resulting  wall 
teoperature  must  he  a  constant.  We  have,  therefore,  demonstrated,  the 
consistency  of  the  assumed  boundary  conditions. 

It  nay  be  of  interest  to  Investigate  the  validity  of  the  modified 
Reynolds  analogy  (which  holds  for  solid*wall  conditions)  under  the  new 
conditions  when  mass  transfer  Is  employed.  If  the  analogy  is  valid,  the 
parameter  (Ojj  Pr®/^)  t  (0^/2)  would  equal  unity.  As  seen  in  Pig.  0, 
considerable  departure  occurs  and  ve  conclude  that,  at  least  for  this 
oonstant'preperty  laminar  flow,  the  modified  Reynolds  analogy  does  not 
apply  In  the  presence  of  mass  transfer. 

AMALYSIS  OP  BARON^^* 

Baron  was  successful  in  Introducing  the  Influence  of  variable 
physical  properties  on  mass-transfer  cooling,  while  at  the  same  time 
beeping  the  n\md>er  of  Independent  parameters  at  a  rilnlnum.  To  acooopllsh 
this,  he  adopted  an  approach  similar  to  that  used  earlier  by  Chapman  and. 
RUbesln^^^^  for  the  solid  flat  plate.  Baron  Introduces  the  so-called 
Obapman-Ruhesln  constant,  C,  for  the  product  of  the  density  and  vlsooslty 
of  the  air  only.  Re  notes: 


where 


CX 


(17) 


C  m  (p«Uft)/(p„u.)  »  Chapman-Rubesln  constant  for  air  only 

<i  G  G  6  rnmmmm 


{Cf/2) 
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Fig.  8  —  Effect  of  moss  transfer 
on  modified  Reynolds  analogy 
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-EiL. 


2  L 


The  factor,  C,  la  Independent  of  the  concentration,  while  the  other 
factor,  "K,  la  a  function  of  both  temperature  (throu^Sh 
concentration.  However,  for  heliisn  and  carbon  dioxide,  the  gases 
ocnsldered  by  Baron,  the  viscoalty  term  n/pg  la  relatively  inaensltive 
to  temperature  and  is  primarily  dependent  on  the  concentration'*^.  Conse¬ 
quently,  Baron  asaumea  that  X  la  a  function  only  of  concentration. 

Using  this  aasuD^tion,  Baron  presents  two  approaches: 

Approach  1 

Baron  introduces 


u  - 

p  ay 


d'i' 


y  u  » 


dj> 

dx 


'I'  »  X  UgC  f  (17) 


(18) 


Using  these  transformations  along  with  the  additional  assumption  that  the 
Schmidt  number  is  a  function  only  of  concentration  (an  Independent  cheek 
shows  this  to  be  a  realistic  assumption),  Baron  obtains  a  set  of  three 
ordinary  differential  equations,  similar  in  form  but  more  con^lex  than 
Eqs.  (12)  -  (16).  The  net  result  of  these  substitutions  is  that  none  of 


*This  assimqjtion  is  valid  for  hydrogen  as  well. 
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the  tenna  appearing  in  the  momantum  and  diffusion  equations  are  dependent 
on  temperatiare  and  consequently  these  tvo  equations  may  be  solved 
simultaneously,  but  Independently  of  the  energy  equatloni  Using  this 
approach  Baron  obtains  the  velocity  profiles,  the  concentration  profiles 
and  the  skin -friction  parameter,  e^  ^u^x/v^C,  as  a  function  only  of  the 
mass-transfer  parameter  (p^v^)/(p^Ug)  .  Any  teinperature 

effect  is  completely  contained  in  the  constant  C. 

The  energy  eqriation  romailns  to  be  solved,  and  Baron  reports  that 
all  the  coefficients  appearing  in  this  equation  were  only  mildly  affected 
by  tenperature,  allowing  the  assunption  that  sdl  coefficients  are 
functions  only  of  concentration.  This  is  a  considerable  slnpliflcatlon 
for  nov  the  energy  equation  is  linear  in  tesperature,  since  all  coeffi¬ 
cients  are  known  functions  of  the  dimensionless  parometer,'*; ,  by  virtue 
of  the  previously  obtained  solutions  of  the  momentum  and  dlff\)slon 
•qmtlons.*  As  In  the  constant-property  situation,  the  resulting 
energy  equation,  being  linear,  nay  be  solved  by  first  treating  the  non- 
disslpatlve  case  and  then  adding  the  adiabatlc-vall  solution.  It 
follows  that  the  low- speed  non-dissipative  heat-transfer  coefficients 
nay  be  used  for  the  hie^-speed  case  if  the  adiabatic -wall  temperature 
replaces  the  free-strean  tenperature  in  the  definition  of  the  heat- 


'M’ There  appears  to  be  an  error  in  Baron's  final  energy  equation  in 

Aef .  2.  This  reevilts  when  he  replaces  rather  than 

rt  e  e  p.;  e 

the  oerreet  ei^gpreeslen  (7.  -  1)  0^  .  Consequently,  the  left  hand 

side  of  Baron's  energy  equation  61X5  should  be  multiplied  through  by 

0  o/c_  to  get  the  correct  form. 

P2'  P- 
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transfer  coefficient 


q  -  h(T^  -  T^)  (19) 

The  recovery  ten^erature,  determined  from  the  adiabatic -wall 
solution.  Is  reported  In  terms  of  a  recovery  factor,  r,  which  for  a 
given  Injected  gas  Is  a  function  only  of  the  mass -transfer  parameter 
/  P-U  )  ^ux/i/ C  .  In  addition  to  the  recovery  factor,  Baron 
also  presents  dimensionless  heat-transfer  coefficients  for  two  binary 
systems,  hellum-alr  and  carbon  dioxide -air  mixtures. 

The  second  approach  used  by  Baron  Is  not  as  realistic  as  the  above 
and  will  be  mentioned  only  briefly.  In  this  case  he  asstimed  that  \  is 
a  constant  to  be  evaluated  at  wall  condition  and  the  following  trans¬ 
formations  are  then  applied  to  Bqs.  (S)  •  (8)s 
Approach  S 


Baron  obtained  seme  representative  solutions  for  this  slurpllf led 
ease  and  ooopared  them  with  the  more  realistic  case  outlined  above.  This 
comparison  la  shown  for  helium  Injection  In  Fig.  9  and  the  agreement  is 
only  fair. 


Heot-tronsfer  parameters,  Cu/C^o  o"*! 


Fig.  9 — Helium  results  of  Baron 
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In  a  later  publication^^ ^  Baron  generalizes  his  analysis  to  include 
the  influence  of  pressure  gradients.  In  addition,  he  presents  a  summary 
of  his  flat-plate  results.* 

ANALYSIS  OF  ECKERT  AND  CO-WORKERS^ 5,  6) 

Before  attempting  the  complete  binary  problem  including  heat  transfer, 
Eckert  and  Schneider  first  solved  the  isothermal  case  with  hydrogen  as 
the  injected  gas.  The  physical  properties  were  allowed  to  vary  with 
concentration  and  the  methods  outlined  in  Ref.  14  were  used  to  calculate 
the  variation.  The  resulting  velocity  distributions  and  skin  friction 
are  con^jared  to  the  constant  property  results  in  Figs.  9  and  10,  This 
con^arlson  reveals  that  unstable  S-shaped  velocity  profiles  occur  at 
relatively  low  values  of  the  dimensionless  mass  transfer  when  compared 
with  the  constant  property  situation.  The  greater  Influence  of  a  light 
gas  on  skin  friction  is  obvious  in  Pig,  11,  where  at  the  same  values  of 
dimensionless  mass  transfer  considerably  lower  skin  friction  occurs  for 
the  hydrogen  injection.  We,  therefore,  conclude  that  the  light  gas  is 
more  effective  in  reducing  skin  friction  but  on  the  other  hand  is  more 
de-stabilizing  to  a  laminar  boundary  layer.  This  will  be  demonstrated 
in  a  later  section. 

In  the  analysis  of  the  binary  system  including  heat  transfer, 
again  using  hydrogen  as  the  maos -transfer  medium,  Eckert  and  his 
colleagues  used  the  following  transformation  in  dealing  with  Eqs.  (2)  -  (8): 


*Care  should  be  taken  in  using  this  reference  since  there  is  some 
confusion  in  nomenclature.  The  mathematical  development  utilizes  a  some¬ 
what  different  transformation  from  that  used  by  Baron  in  Ref.  2,  although 
the  figures  are  all  shown  in  terms  of  the  original  variables  of  Ref.  2. 
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The  transformed  equations  were  then  solved  using  the  best  available 
property  Information  for  hydrogen-air  mixtures >  These  exact  solutions | 
which  were  obtained  using  an  iterative  procedure  on  an  ERA.  1103  elertronie 
cotqputerj  are  valid  only  for  the  specific  conditions  selected*  (See 
rigs*  2  and  3*)  An  exaople  of  these  results  is  shown  in  Fig.  12  where 
the  dimensionless  heat-transfer  coefficients  are  shown  for  eero  Mach 
number  for  two  different  wall-tesperature  conditions  with  the  free 
stream  at  392°R,  Additionally^  for  Mach  12,  the  Nuseelt  number  and 
recovery  teiqperature,  both  of  which  must  be  known  to  determine  the  heat 
transfer,  are  shown  for  a  set  of  specific  conditions*  In  every  case  we 
find  a  considerable  reduction  in  the  heat  transfer  when  only  small 
amounts  of  hydrogen  are  transferred  away  from  the  wall  into  the  boundary 
layer*  The  effectiveness  of  lied^t-gas  injection  in  decreasing  the  heat 
trcmsfer  in  a  high-speed  boundary  layer  is  obvious  from  this  flgicre* 

ANALYSIS  OF  SZIKLA3  AMD  BAMAS^*^^ 

Sslklas  and  Banas  report  solutions  for  a  number  of  different 
coolants:  hydrogen,  helim,  water  vapor,  and  air*  In  arriving  at  the 
final  form  of  the  energy  equation  they  assumed  that  the  specific  heat 


0.25 


Effect  of  hydrogen  injection  on  laminar  heat  transfer 
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ratio,  7,  is  the  same  for  the  injected  coolant  as  for  the  main  stream  gas 
(7^  =  7g).  They  then  use  the  standard  Blasivis  -  transformation 
essentially  as  given  in  Kq.  (ll).  Iji  obtaining  their  solutions,  the 
physical  properties  (inclviding  specific  heat)  were  allowed  to  vaxy  with 
both  ten^erature  and  concentration.  Methods  of  kinetic  theory  were  used 
in  the  determination  of  these  properties.  As  was  the  case  with  the 
results  of  Eckert,  the  results  of  Szlklas  and  Banas  are  applicable  only 
to  the  specific  conditions  imposed  in  the  analysis.  Representative 
results  for  the  helium  study  are  given  in  Fig.  I3,  where,  again,  large 
reductions  in  heat  transfer  accompany  small  mass -transfer  rates. 

ANALYSIS  OF  QROSS^^^ 

The  isothermal  laminar  binary  boundary  layer  on  a  flat  plate  was 

(8) 

investigated  by  Gross'  ^  for  three  different  Injactant  gases:  hydrogen, 
carbon  dioxide,  and  iodine  vapor.  To  obtain  a  solution  of  the  governing 
equations  (Eqs.  (2),  (3),  and  (3))4  Gross  used  the  standard  Blaslus 
transformation  to  arrive  at  a  system  of  ordinary  differential  equations. 
These  were  then  solved  using  a  Runge-Kutta  numerical  method  \riLth  the 
aid  of  a  high-speed  electronic  cosqiuter. 

The  resulting  values  of  the  skin-friction  ceefflclent  for  the  three 
gases  investigated  are  shown  in  Fig.  l4.  These  results  demonstrate  that 
the  addition  of  a  heavy  gas  such  as  iodine  vapor  (molecular  weight  253.0) 
is  much  less  effective  in  reducing  the  skin-friction  coefficient  than  the 
addition  of  the  lighter  gases. 


Dimensionless  moss  fronsfer, 


vRe7 


Vw 
P*  Lie 


3  —  Variation  of  laminar  heat  tranfer 
with  helium  injection 


Dimensionless  moss  tronsfer, 

re  ^9  ^ 

Fig,  14  —  Friction  coefficient  os  a  function  of  the  biowing  parameter  (8) 
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IV.  GEWTOAIJZRD  PRKOKNTATION  OF  LAMINAR  FLAT  PIATK 
BINARY  BOUNDARY  LAYER  RESULTS 

It  is  our  goal  In  this  section  to  develop  a  generalized  presentation 
which  may  be  conveniently  utilized  by  design  engineers  for  predicting 
skin  friction  and  heat  transfer  in  the  presence  of  mass>transfer  cooling 
for  laminar  flow  over  surfaces  with  zero  pressure  gradient.  The  avail¬ 
able  analytical  solutions  briefly  described  above  are  used  as  the  basis 
for  the  generalization.  Since  heat  transfer  and  skin  friction  for  solid 
sux^aces  in  the  absence  of  mass-transfer  cooling  can  be  calculated  at 
the  present  time  with  a  measure  of  confidence,  the  approach  adopted  here 
is  to  present  the  correction  factors  which  must  be  applied  to  such  solid 
wall  calculations  to  account  for  the  effect  of  mass  addition.  Thus  the 
normalized  skln-frlctlon  coefficient  and  heat  transfer  will  be  given  as 
Oj/ejo>  80^  where  the  subscript  zero  Implies  that  the 

q.uantlty  Is  to  be  evaluated  for  the  same  free-stream  and  wall-ten^>erature 
conditions,  neglecting  the  Influence  of  mass  transfer. 

It  was  found  that  these  normalized  results  for  any  one  gas  could  be 
presented  as  a  unique  function  of  the  mass -transfer  parameter  proposed 
by  Baron,  (p^v ^p^u^ )  ^u^  x/  ,  provided  that  the  Chapman-Rubesln 
constant,  0,  was  evaluated  at  the  so-called  reference  temperature,  T*, 
given  by  Eq.  (9).  The  success  of  this  generalized  presentation  for  skin 
friction  Is  demonstrated  In  Figs.  15  and  I6  which  apply  to  hydrogen -edr 
and  hellum-alr  binary  systems,  respectively.  It  appears  that  this 
representation  Is  valid  over  a  wide  range  of  wall-temperatiu'e  conditions 
and  free-stream  Mach  numbers.  Tills  conclusion  is  true  for  the  other 
binary  systems  as  well,  and  the  reader  Is  referred  to  Appendix  C  for 
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Fig.  16  —  Effect  of  moss  transfer  on  skin  friction 

Helium— air 


RM-2516 

40 


verification  of  this  conclusion.  A  summary  of  available  skin-friction 
results  for  six  gases  is  given  in  Fig.  I7. 

This  same  presentation  was  successful  in  correlating  the  normalized 
Stanton  number,  with  the  final  results  ns  shown  In  Pig.  18.  Since  the 
determination  of  the  actval  heat  transfer  requires  the  knowledge  of  the 
recovery  factor  as  well  as  the  Stanton  number,  the  normalized  recovery 
factor  is  shown  in  Fig.  19.  It  may  be  seen  that  some  disagreement 
exists  for  the  light-weight  gases.  Apparently  the  recovery  factor  is 
somewhat  more  sensitive  to  the  different  assumptions,  particularly 
physical  property  variations,  adopted  by  the  various  Investigators.  The 
effect  on  the  final  heat-transfer  prediction  of  this  disagreement  in 
recovery  factor  is  reduced  if  the  normalized  heat  transfer  is  directly 
considered  rather  than  the  Stanton  number.  The  resulting  normalized 
heat  transfer  is  shorn  in  Figs.  20  through  22.  It  should  be  pointed  out 
that  some  effect  of  the  disagreement  in  recovery  factor  is  still  present 
in  this  presentation  through  the  presence  of  O’*'.  However,  for  practical 
applleatlons  the  wall  tenperature  in  general  will  be  markedly  lower  than 
the  recovery  tenperature}  for  this  situation,  inspection  of  the  defining 
equation  (9)  for  T*  reveals  that  any  uncertainty  in  the  recovery  tem¬ 
perature  has  only  a  minor  effect  on  the  reference  tenperature  Itself, 
and  consequently,  only  a  minor  effect  on  C’*'. 

It  is  apparent  from  Figs.  17  and  22  that  the  light  gases  are  much 
more  effective  than  the  heavier  gases  in  reducing  heat  transfer  and  skin 
friction.  Inspection  of  these  figures  indicates  that  the  normalized 
skin-friction  coefficient  and  heat  ti’ansfer  q/q^  vary  linearly 

with  the  dlmennlonless  mass-transfer  parameter  for  all  the  gases  shown. 
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In  particular,  the  results  for  air  into  air  may  be  expressed  by  the 
following  two  equations: 

q/q^,  -  1  -  1.82  {(PyVyp^Ug)  CyRe^)]- 

(22) 

®f/®fo  - 1  ‘  2*08  {(p^vyp^u^)  (yReyc*)} 

It  is  of  tec^lcal  Isqiortanoe  to  determine  whether  the  ether  gases 
oan  be  made  to  agree  with  these  equations  by  the  slnple  expedient  of 
multiplying  the  Baron  dimensionless  mass-transfer  parameter  by  a  molecular 
weight  ratio  (mg/m^),  raised  to  a  constant  exponent.  This  question  can 
be  answered  by  plotting  the  dimensionless  mass  transfer  versus  the  mol¬ 
ecular  weight  at  a  constant  value  of  q/q^  (or  Of/cjo)*  I® 

aceoBgtllshed  In  Fig.  23  where  It  Is  seen  that  l/3  represents  a  fair 
conqwomlse  for  the  value  of  the  exponent  although  the  very  light  gases 
as  well  as  the  heavier  gases  such  as  Iodine  show  a  significant  departure. 
Recognizing  that  such  discrepancies  do  occur  for  the  heavier  and  light 
gases,  It  nevertheless  appears  that  the  following  equations  represent 
the  heat  transfer  and  skin  friction  reasonably  veil. 

q/q^  -  1  -  1.82-[(m2/m^)^/^  ^Vw/Pe^e^  (<v/^^7^)} 

(23) 

Of/Cfo  -  1  -  2.08  {(mg/m^)^'^  (p^v^p^u^)  (/Reyc*)} 

It  Is  recommended  that  these  equations  be  used  to  predict  the  heat- 
transfer  and  Bkin-frlctlon  performance  for  mass -transfer  cooling  In  a 
binary  laminar  boundary  layer  on  a  flat  plate. 


versus  molecular  weight  of  coolant  gas 
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V.  EFFECT  OF  PRESSURE  GRADIENT  ON  LAMINAR 
MASS -TRANSFER  COOLING 


to  the  present  time,  little  effort  has  been  directed  to  the 

solution  of  the  binary  laminar  boxmdary  layer  equations  vlth  finite 

pressure  gradients.  However i  a  measure  of  the  ixiflusnce  of  a  pressure 

gradient  on  mass-transfer  cooling  can  be  obtained  by  returning  to  the 

constant-property  boundary  layer  model  with  normal  Injection  (air-lnto- 

alr)  since  solutions  have  been  reported  in  this  case  for  wedge-type 

pressure  gradients  (l.e.,  the  free-stream  velocity  Is  described  by 

u_  "Ax®).  Examination  of  these  solutlons^^^^  reveals  that  the 
o 

presence  of  a  favorable  pressure  gradient  leads  to  more  stable  velocity 
profiles  and,  consequently,  It  appears  that  larger  values  of  the  dimen¬ 
sionless  mass-transfer  parameter,  (p^y^p^u^  may  Ts«  obtained 

without  causing  transition  to  turbxilence.  Furthermore,  the  skln-frletlon 
coefficient  remains  finite  in  a  favorable  pressure  gradient,  with  no 
apparent  failure  of  the  boundary  layer  equations  even  for  very  large 
mass-transfer  rates.  However,  the  heat  transfer  does  decrease  to  a 
diminishing  value,  with  the  thermal  boundary  layer  being  displaced  from 
the  wall  surface  toward  the  free  stream.  An  example  of  this  skln- 
frlctlon  and  heat-transfer  behavior  Is  given  In  Fig.  24  for  plane 
stagnation  flow  (  m  l),  and  it  may  be  seen  that  the  heat  transfer 
goes  to  zero  at  a  value  of  the  mass -transfer  parameter  of  approximately 
2  {as  contrasted  to  O.619  for  the  flat  plate),  while  the  skin  friction 
is  still  approximately  UO  per  cent  of  its  eolld-Jwall  value. 

A  direct  conparison  of  the  reduced  heat  flow,  (l/q^j,  for  four 
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different  pressure  gradients  ranging  from  the  zero  pressure  gradient 
flat  plate  to  the  plane  stagnation  flow  is  given  in  Fig.  25.  If  a  com¬ 
parison  is  made  at  a  fixed  value  of  the  dimensionless  mass-transfer 
parameter,  it  is  apparent  that  the  greatest  reduction  in  heat  transfer 
occurs  for  the  zero  pressure  gradient  flat  plate  with  the  least  reduction 
accompanying  the  plane  stagnation  flow.  Therefore,  to  obtain  a  given 
reduction  in  heat  flow,  it  is  necessary  to  go  to  higher  values  of 

the  dimensionless  mass  transfer  as  the  pressure  gradient  increases. 

Since  there  exists  considerable  Interest  in  the  three-dimensional 
stagnation  flow,  Fig.  25  also  presents  the  reduced  heat  transfer  q/q^ 
for  air  injection  into  such  a  region^^*^^. 

It  finally  remains  to  determine  whether  the  relative  position  of 
the  varlo\is  coolant  gases  found  for  the  flat  plate  geometry  is  markedly 
influenced  by  the  presence  of  a  pressure  gradient.  A  recent  analysis 

(it) 

by  hayday'  '  for  hydrogen  injection  into  a  plane  stagnation  flow  leads 
to  the  results  shown  in  Fig.  26  and  close  inspection  suggests  that  the 
molecular  weight  parameter  found  for  the  flat  plate,  i.e.  (mg/nij^)^/^,  is 
approximately  valid  for  the  plane  stagnation  flow. 

As  a  result,  it  is  suggested  that  alr-into-alr  resiolts  be  used  to 
predict  the  effect  of  preos\are  gradient  on  laminar  mass-transfer  cooling; 
if  a  coolant  other  than  air  is  transferred  into  the  boundary  layer  the 
relative  effectiveness  of  the  various  coolants  is  to  be  estimated  fran 
the  flat-plate  results.  This  is  siaqjly  accomplished  by  ^^sing  Fig.  25, 


changing  the  abscissa  to  read 


(PwV^e^e) 


1.2 


Constant  physical  properties 
Pr  =  0,7 


Fig.  25 — Effect  of  pressure  gradient 
on  mass -transfer  cooling 
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VI.  STABILrn:  of  the  laminar  binary  boundary  layer 

THEORETICAL  8TUDIE3 

The  Tollmlen-Sehliehtlng  theory  of  email  dleturbanoeo  hae  been 
widely  used  to  Investigate  the  stability  characteristic's  of  laminar 
boundary  layers  on  solid  surfaces  and  recently  has  been  applied  to 
binary  boundary  layers.  Essentially^  this  theory  Is  a  perturbation 
analysis  of  the  Navier-Stokes  equations.  Secular  relationships  between 
the  coefficients  of  the  linearized  perturbation  equations  define  a 
region  In  which  disturbances  of  a  given  wave  length  are  damped  or 
excited.  A  damped  disturbance  suggests  that  eventually  turbulence  will 
occur. 

Lin  eind  Lees^^^^  have  presented  a  stability  analysis  for  a  single- 
component,  compressible,  bovmdary  layer  flow  over  a  flat  plate.  The 
computational  effort  required  to  define  the  coefficient  relationships 
is  rather  foraldable.  Consequently,  some  approximations  have  been 
suggested  by  Lln^^^^  for  the  Incompressible  flow  case  and  by  Lees^^^^ 
for  the  compressible  flow  case  which  will  be  adequate  to  predict  the 
minimum  critical  Reynolds  number  below  which  the  Tollmlen-Schlicbting 
dlsturbancee  of  all  wave  lengths  are  damped;  hence,  the  flow  should  be 
stable  for  Reynolds  numbers  below  this  value. 

The  simplified  Lin  approach  yields  the  following  criterion  for 
determining  the  mlnlmxm  critical  Reynolds  number; 

^*“^crj  "  d(y/8*)  '^i  (24) 


where 


6* 
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6 


*  o 


dy 


where  la  determined  by  examining  the  apeelflc  laminar  velocity 
profile  which  appllea  to  the  condltlona  under  Investigational  is 
that  locaJL  velocity  existing  at  the  position,  y,  such  that  the  following 
equation  Is  satisfied] 


0.58  -tK^) 


■2(du/dy)^j^ 

'u(y^)  (d®u/dy®)y^" 

■  u(y^)  “  5 

(du/cly)^y 

(25) 


Recently  Gross  has  suggested  that  the  Lin  equations  (2h)  and  (S^) 
be  utilized  as  a  first  approximation  for  predicting  instability  in  a 
binary  laminar  boundary  layer. A  comparison  of  the  predicted  minimum 
Reynolds  number  for  relative  destabilizing  effect  should  Indicate 
direction  Insofar  as  the  different  injection  materials  are  concerned. 

In  the  application  of  Lin's  single  con^onent  approximation  to  a  binary 
system,  an  Isothermal  boundary  layer  was  assumed*  Furthermore,  the 
amount  of  the  Injected  material  was  considered  to  be  so  small  that  (a) 
the  boundary  layer  assunptlons  were  not  Invalidated,  emd  (b)  the  effects 
of  the  injected  material  would  influence  the  solution  only  In  the 
variation  of  the  property  parameters.  Essentially,  this  means  that  the 
distxsrbonoe  equatlcnu  for  the  velocity  and  concentration  are  Independent. 
The  velocity  and  density  profiles  found  by  solving  the  momentum  and 
diffusion  equations  were  substituted  into  Eqs.  (S4)  and  (2^)  yielding 
the  results  shown  In  Fig,  27.  It  la  apparent  that  In  an  Isothermal 
binary  boundary  layer,  injection  of  any  material  will  be  destabilizing. 
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but  this  effect  tends  to  decrease  as  the  molecular  weifjht  of  the 
material  Increases.  The  hydrocen,  then,  appears  to  be  more  destabnizlng 
than  the  heavier  gases,  carbon  dioxide  and  Iodine.  This  la  the  result 
primarily  of  the  higher  velocity  necessary  to  supply  the  low  density 
gas  in  sufficient  amounts  to  satisfy  the  blowing  parameter  relationship. 
The  Injection  of  hydrogen  at  the  wall  is  somewhat  analagous  to  heat 
transfer  from  the  wall  In  that  It  redtices  the  density  in  the  neighborhood 
of  the  wall.  This  density  reduction  encoiarages  a  corresponding  decrease 
In  the  skln-fr lotion  coefficient.  It  has  been  shown  in  earlier  sections 
that  hydrogen  la  the  moat  effective  agent  for  mass -transfer  reduction  of 
heat  transfer  and  akin  friction.  It  Is,  however,  also  the  most  desta¬ 
bilizing.  Iodine  vapor  and  carbon  dioxide,  which  are  less  efficient  in 
reducing  heat  transfer  and  skin  friction,  show  much  less  destabilizing 
action.  As  the  first  approximation,  Pig.  27  may  be  used  In  conjunction 
with  Fig.  22  to  determine  the  optimum  injectant  to  effect  a  given  wall 
temperature  condition  while  maintaining  a  maximum  stability.  A  study  is 
now  under  way  to  obtain  some  very  accurate  compressible  multi -conponent 
boundary  layer  Information  which  will  be  acceptable  for  a  more  rigorous 
stability  analysis. 

The  more  complex  problem  of  the  compressible,  binary,  boundary  layer 

(21) 

flow  has  been  treated  by  E.  E.  Covert  using  a  Tollmlen-Schlichting 
perturbation  analysis.  He  showed  that  the  perturbation  equations  may 
be  solved  in  exactly  the  same  manner  as  for  a  single  component  gas 
except  for  the  implicit  effects  of  the  foreign  gas  on  the  local  flow 

(0 ) 

properties, thereby  Justifying  the  approach  of  Gross.'  '  This  procedure 
holds  for  relatively  small  injection  rates,  where  the  boundary  layer 
assumptions  are  still  valid.  Covert  found  that  the  stability 
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charncterlstlcB  could  be  oalevilated  for  an  n~cemponent  boundary  layer 
using  the  steady-state  velocity  and  density  profiles.  For  a  given  range 
of  Mach  number,  one  may  establish  a  boundary  laj^r  completely  stable  to 
two-dimensional  subsonic  disturbances  sin^ily  by  cooling  the  surface. 
Covert  calculated  values  of  the  iraill  ten^erature  required  for  con^lete 
stability  over  a  range  of  Mach  ntanbers  with  the  mass -transfer  rate  as 
another  parameter  (see  Fig.  28).  It  should  be  noted  that  his  calcu¬ 
lations  are  based  on  the  solution  of  the  Invlscid  disturbance  equations. 
Helium,  carbon  dioxide,  and  air  injection  were  considered. 

In  general,  it  may  be  noted  for  all  three  Injectants  that  an 
increase  in  the  injection  rate  normally  results  in  lower  sxirface 
temperatures  at  the  re-entry  speeds  of  present  day  Interest.  This 
decrease  in  the  surface  teirperature  tends  to  stabilise  the  laminar 
’•ovmdary  layer.  This  should  he  particularly  true  for  the  injection  of 
huavy  molecviles  since  these  will  result  in  the  same  effect  as  cooling 
on  a  solid  wall,  namely  increasing  the  density  near  the  surface.  How¬ 
ever,  it  is  obvious  that  the  act  of  injection  Introduces  a  disturbance 
in  the  boundary  layer  by  decelerating  the  flow.  This  will  certainly  be 
destabilizing  at  high  injection  rates.  Therefore,  we  may  expect  a 
counterbalancing  of  stabilizing  and  destabilizing  effects.  The  net 
result  will  depend  upon  the  efficiency  of  the  injected  material  in 
cooling  the  boundaiy  layer.  Figure  29  shows  the  minimum  tenperature 
ratio  for  the  injection  of  helium  as  a  function  of  Mach  number  with  the 
dimensionless  mass  transfer  rate  as  a  parameter  as  determined  by  Covert. 
It  Is  Indeed  so  that  for  a  Mach  number  below  3.5,  the  injection  of  this 
light  gas  reduces  the  stability  character  of  a  single  conponent  layer. 


Fig.  28 — Effect  of  helium  injection  on  cooling  reguired  for  complete  boundary  layer  stobllit 
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Fig.  29 — Effect  of  mass  transfer  on 
boundary  layer  transition 
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However^  as  the  Mach  number  increases,  the  adiabatic  vail  temperature 
rises  very  rapidly  causing  a  high  enthalpy  region  near  the  vail.  The 
helium  with  Its  high  specific  heat  injected  into  the  hot  boundary  layer 
acta  as  a  coolant  decreasing  the  enthalpy  near  the  wall,  and  the  over¬ 
all  effect  is  to  Increase  the  minimum  wall  ten^iorature  ratio.  An 
interesting  conclusion  of  Covert  is  that  no  amount  of  air  or  carbon 
dioxide  injection  will  Increase  the  minimum  temperature  ratio.  This 
agrees  with  the  result  shown  previously  by  the  isothermal  Reynolds 
number  calculations. 

EXPERIMENTAL  STUDIEa 

(7) 

Sziklas  and  Banas' ' '  have  reported  results  of  experimental  work 
with  a  T-inch  17* deg  apex  angle  porous  cone  operating  at  Mach  2.6.  Dae 
injected  coolants  were  helium  and  air.  Velocity  profiles  at  various 
injection  rates  are  given.  The  helium  data  for  (py* V^’e'^e ^ '^^x  ^ 
indicated  profiles  which  agreed  with  theoretical  calculations,  but  at 
a  value  of  0.14  this  agreement  was  no  longer  true  and  the  profile  was 
transitional.  The  corresponding  case  for  air  showed  agreement  tip  to 
(p^v^PgUg)  »  0.02,  but  the  last  injection  rate  again  gave  a 

transitional  profile.  These  results  indicate  that  mass  injection  rates 
of  sufficient  magnitude  to  produce  significant  cooling  effects  do  not 

show  a  pronounced  tendency  toward  turbulence. 

(2'5) 

Leaden,  Scott,  and  Anderson  '  used  a  20-deg  angle  cone  at  a  Mach 
number  of  5«  The  Injcctants  were  helium  and  air.  The  non-ideal 
permeability  (and  consequent  use  of  a  modified  blowing  parameter)  of 
the  model  and  large  scatter  in  the  data  make  interpretation  difficult. 
Nevertheless,  it  was  shown  that  the  heat-transfer  results  of  the  porous 
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model  (zero  blowing)  were  quite  similar  to  those  for  a  smooth  model. 

The  transition  Reynolds  number  with  zero  blowing  (  ~  2.5  x  10^)  on  the 

porous  model  was  about  half  as  Isirge  as  that  for  a  similar  smooth 

model.  Finally,  the  authors  point  out  that  the  boundary  layer  was 

quite  stable  to  air  and  helium  injection. 

(24) 

Scott,  Anderson,  and  Elgin  have  made  measurements  on  the  effects 

of  helium  and  air  injection  on  boundary  layer  transition  on  a  l6-deg 

porous  cone  modal  at  a  Mach  number  of  3.7.  They  determined  transition 

position  with  the  use  of  shadowgraph  photographs.  The  dimensionless 

transition  Reynolds  number  was  plotted  against  the  corresponding  value 

of  the  injection  parameter,  as  shown  in  Fig.  29.  The  percentage  decrease 

of  Re^/Re^^  is  not  large  for  Injection  rates  which  produce  significant 

cooling  effects.  In  addition,  the  order  of  the  mass-transfer  rates  and 

the  rate  of  decrease  between  helium  and  air  is  similar  to  the  Isothermal 

(8) 

theoretical  results  obtained  earlier.  ' 
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VII.  MABS-TRAMSFER  COOLINO  IN  A  TURBULENT  BOUNDARY  LKCER 

thllke  laaloar  flow,  it  la  not  possible  to  solve  rigorously  the 
governing  equations  for  tuijbulent  binary  boundary  layer  flow  in  a 
direct  nanner.  Rather,  a  variety  of  indirect  approaches  have  been 
suggested  ranging  from  the  slnqple  film  theory  to  more  complex  mixing 
length  analyses.  The  simple  film  theory  or  Couette  flow  model 
conceptually  replaces  the  actual  boundary  layer  by  a  film  of  constant 
thloknesB,  neglecting  any  variations  in  the  flow  direction  as  eoutpared 
with  those  across  the  fUm.^^^'  Although  this  approach  does  not 
yield  all  details  of  the  flow,  it  does  predict  the  over-all  effect  of 
mass  transfer  on  skin  friction  and  heat  transfer  quite  well.  More 
detailed  mixing  length  analyses  which  require  engtlrical  data  to 
aeoonpllsh  a  final  solution  have  been  advanced  for  air  injection  by 
Dorranee  and  Dore,^®"^^  Rubesln, and  van  Driest. The  first  two 
references  era  applieahle  to  the  flat  plate  geometry,  with  Dorranee  and 
Dore  raatrloting  their  analysia  to  a  Prandtl  number  of  unity  and 
aaaumlng  the  turbulent  boundary  layer  extends  down  to  the  surface. 

The  analyses  of  Rubesln  apply  to  a  Prandtl  number  of  0.7^  consider 
viscous  dissipation,  and  allow  for  a  laminar  sublayer  and  a  turbulent 
core.  The  incompressible  analysis  of  van  Driest  considers  the  laminar 
sublayer  and  turbulent  core  using  somewhat  different  assumptions  than 
those  of  Rubesln  for  the  flat-plate  geometry,  but  in  addition,  van 
Driest  extends  the  analysis  to  apply  to  air  injection  into  a  turbulent 
boundary  layer  in  the  region  of  either  a  two-dimensional  or  three- 
dimensional  stagnation  point.  Tlie  results  of  Rubesln  and  van  Driest 
are  essentially  in  agreement  with  each  other  for  the  flat-plate 
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geometry  and  are  somewhat  higher  than  the  prediction  of  Dorranoe  and 
Dore,* 

A  ntimber  of  experimental  restilte  are  reported  for  air  injection 
into  a  turhtilent  boundary  layer  on  zero  preaeuro  gradient  surfaces,'  ' 
31-55)  thereby  providing  a  check  on  the  above  predictions.  Figure 
30  CQn^ares  the  available  heat-transfer  results  In  tbs  form  of  Stanton 
numbers  with  the  several  analyses.  It  may  be  seen  that  the  three 
analyses  do  not  differ  very  markedly  although  Dorrance  and  Dora  are 
somewhat  lower.  The  experimental  data  covering  a  range  of  Mach  nuobers 
from  0  to  3  show  considerable  scatter  but  nevertheless  the  trend  is 
correctly  predicted  by  all  three  analyses.  Similar  results  for  skin 
frietlon  are  shown  in  Fig.  31  and  the  bulk  of  the  data  tend  to  agree 
with  the  prediction  of  Rubesln  and  van  Driest.  The  only  conflicting 
evidence  is  the  experimental  data  of  Mlckley  and  Davis' which  show 
skln-frintlon  values  much  lower  than  the  other  Investigators  —  even 
lower  than  predicted  by  the  film  theory.  Since  the  majority  of  the 
experimental  data  favor  Rubesln  and  van  Dr lest >  it  appears  at  the 
present  time  that  this  represents  the  more  reallstie  and  certainly  the 
more  conservative  estimate  of  skin  friction  and  beat  transfer  in  the 
presence  of  air  injection. 

The  effect  of  pressure  gradient  on  heat  transfer  and  skin  friction 
in  the  presence  of  mass  transfer  in  turbulent  flow  is  apparently 
minimized  when  the  results  are  presented  in  the  coordinates  shown  in 
Figs.  30  and  31.  ansQysls  of  van  Driest  indicates  that  the  same 


♦Dorrance  discusses  this  discrepancy  Vhich  arises  from  the  form  he 
assumes  for  his  velocity  distribution. '30) 


Effect  of  moss  transfer  on  heat  transfer,  turbulent  flow 


Predicted  and  measured  skin  friction  for  air  injection  into 
a  turbulent  boundary  layer  on  a  fiat  plate 
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curve  shown  for  the  flat  plate  may  also  be  used  for  turbulent  flow  In 
the  region  of  two-  and  three-dimensional  stagnation  points. 

The  calculation  of  heat  transfer  req^llre8  the  knowledge  of  the 
recovery  factor  in  addition  to  the  above-cited  Stsmton  number.  The 
theory  of  Rubesin  predicts  that  there  is  no  effect  of  air  injection  on 
the  recovery  factor  for  turbiHent  flow.  However,  this  is  in  dlfjagree- 
ment  with  experimental  evidence. Since  seme  of  the 
assusqptlona  used  in  the  analysis  of  Rubesin  have  not  been  confirmed, 
greater  weight  must  be  given  to  the  experimental  results.  These  are 
summarized  in  Fig.  32  and  the  ouznre  drawn  through  the  data  may  be  used 
as  a  guide  in  predicting  the  recoveiry  factm*  for  air  injection  into  a 
turbi..  ut  boundary  layer. 

For  foreign  gas  addition  to  a  turbulent  boundary  layer  on  a  flat 
plate,  the  skln-frletion  results  of  Pappas''^  '  and  Ward  and  Bhmon* 
are  available  and  are  represented  in  Fig.  33,  These  predictions,  valid 
for  low  Mach  numbers,  demonstrate  the  advantage  of  the  lightweight 
gasee,  helium  and  hydrogen,  in  reducing  ekln  friction.  The  heavier 
gases,  freon  and  teflon,  are  obviously  less  effective  in  reducing  skin 
friction,  but  eturprlslngly  their  performance  is  not  appreciably  poorer 
than  that  of  air. 

For  heat  transfer,  the  mixing  length  analysis  of  Rubesin  and 
(37) 

Pappas  predicts  the  performance  shown  in  Fig.  34  for  light-gas 
injection  Into  a  turbulent  air  boxuidary  layer  on  a  flat  plate.  The 
limited  available  experimental  data  for  helium-air suggest  an 
even  greater  reduction  than  predicted  the  analysis.  The  predicted 
recovery  factor  shows  no  change  from  the  solid  wall  value  when  light 


*T,  E,  Ward  and  D.  B.  Harmon,  personal  communication,  March  1959 


Fig.  32  —  Recovery  factor  for  air  injection  into  turbulent  boundary  layer  on  a  flat  plate  (31,33) 


Rg.  34  —  Effect  of  moss  transfer  on  heat  transfer,  turbulent  flow 


BM-2916 

73 


gases  are  Injected,  into  a  turbulent  boundary  layer,  and  at  the  present 
tine  the  e^gperlmental  data  are  Insufficient  to  justify  any  definite 
ooneluslon. 

To  detezmlne  the  Influence  of  moleeviLar  welc^t  of  the  coolant  gae 

on  the  heat  transfer  and  skin  friction,  the  mass-transfer  parameter  vas 

plotted  against  the  molecular  velght  at  a  constant  value  of  the  reduced 

paramster,  results  are  shown  in  Fig.  3^ 

0 

where  It  Is  seen  that  the  molecular  wel^t  appears  to  play  a  more 
sensitive  role  for  ...j  light  gas,  while  the  heavier  gases  are  less 
sensitive  to  this  property.  An  estimate  of  the  reduced  skln-frlctlon 
coefficient  and  Stanton  numbers  for  other  coolants  may  be  obtained  from 
the  presentation  In  Fig.  35 • 


Turbulent  mass-transfer  parameter  versus  molecular  vreight  of  coolant  gas 
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VIII.  SUBLIMATION  COOLINQ 

Sublimation  Is  defined  as  the  physical  change  encountered  by  a  Bub> 
stance  In  passing  from  the  solid  to  the  gaseous  phase.  It  Is  characterised 
by  the  absence  of  any  Intermediate  liquid  phase.  Reference  to  a  pressure* 
ten^erature  phase  diagram  makes  it  apparent  that  sublimation  processes 
occur  at  relatively  low  presstires  and  temperatures.  The  tenqterature  at 
vhich  the  vapor  pressure  of  the  solid  equals  the  total  pressure  of  the 
gas  phase  in  contact  with  it  is  defined 
as  the  sublimation  point.  The  snow  point 
is  defined  as  that  temperature  at  vhich 

P 

the  vapor  pressure  of  the  solid  is  equal 
to  the  partial  pressure  of  the  substance 
in  the  gaseous  phase.  It  is  clear  then 

that  our  interest  lies  In  the  snov  point  and  not  in  the  Bubllnatlon 
point.  A  foreign  gas  concentration  of  unity  at  the  surface  (l.e.>  sur* 
face  at  the  sublimation  point)  vo\ild  Imply  a  blowing  rate  high  enough  to 
blow  the  boundary  layer  from  the  surface.  In  general,  the  vapor  pressure 
of  the  solid  material  will  be  a  function  only  of  its  temperature  for 
equilibrium  conditions. 

The  BUbllmatlon-coollng  system  has  already  been  described  briefly 
in  the  Introduction.  The  phase  chEinge  occurring  at  the  vail  in^oses  an 
additional  restraint  on  the  boundary  layer  equations  in  the  form  of 
another  boundary  condition  vhich  relates  the  wall  ten^erature  to  the 
concentration  of  foreign  material  by  the  Clausius -Clapeyron  equation 
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- ; -  -  ~  exp  (  -  M  /RT  )  (26) 

m  -  (m  -  1)  Pe  ® 

Y  m  concentration  of  sublimating  material  at  the  vail 

V 

m  • 

B  m  constant  of  integration 

p  ■  pressure  at  the  edge  of  boundary  layer 

■  heat  of  sublimation  per  mole 
s 

R  ■  universal  gas  constant 

The  introduction  of  this  boundary  condition  lollies  that  a  set  of  unique 
values  of  the  boundary  layer  functions  exist  once  the  external  pressure, 
teQg>erature,  Mach  number,  and  the  stirfaoe  material  have  been  specified. 
Changing  the  external  conditions  then  automatically  defines  a  new  set  of 
functions.  Consequently  ve  may  regard  the  subliming  systenw  as  self* 
regulating  since  the  conditions  at  the  vail  respond  automatically  and 
viniquely  to  changes  in  free-stream  conditions.  The  variation  of  q/q^ 
with  the  mass •transfer  parameter  for  a  COg-alr  mixture  is  shown  in 
Fig.  22.  If  ve  specify  that  solid  carbon  dioxide  is  subliming  into  a 
laminar  boundary  layer  with  fixed,  external  flow  conditions,  then  only 
one  point  on  the  curve  will  satisfy  all  the  boundary  conditions .  Since 
Baron's  parameter  correlates  the  Mach  number  and  wall-temperature  effects 
on  the  heat  transfer,  it  is  possible  to  approxlmat”’  the  variation  of  the 
normalized  heat  transfer,  with  the  mass -transfer  parameter  by  a 

straight-line  relationship.  This  fact  coupled  with  linear  approximation 
to  the  mass-transfer  parameter -wall  concentration  (see  Appendix  B) 
relation  and  the  Clausius -Clapeyron  equation  will  generate  a  simple 
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algebraic  repreaentation  between  the  vmll  temperature  of  a  aubliming 
material  and  the  external  conditions. 

This  equation  was  derived  and  applied  to  systems  In  which  ice, 
carbon  dioxide,  and  ferrous  chloride  are  subliming.  The  results  are 
shown  in  Fig.  36.  The  wall  tenperature  is  quite  insensitive  to  Mach 
number.  On  the  other  hand,  variation  In  altitude  (external  pressixre) 
brings  about  marked  changes  in  the  tenperature  at  the  surface.  TIae 
external  pressure  Is  In  the  denominator  of  the  Claus lus-Clapyron  equation 
and  Is  quite  Inportant  in  determining  the  level  of  the  wall  temperature. 
The  increased  heat  load  brou({ht  about  by  larger  Mach  numbers  Is  absorbed 
by  the  greater  mass-transfer  parameter  (increased  surface  sublimation) 
rather  than  by  large  temperature  increases.  This  accounts  for  the 
insensitivity  to  Mach  number  variation. 

Our  equations  are,  of  eoucrse,  dependent  on  the  condition  of  equil¬ 
ibrium.  As  the  external  pressure  decreases,  the  sublimation  pressure  is 
reached  very  quickly  and  at  extremely  low  wall  concentrations.  This 
inplles  that  only  very  high  blowing  rates  are  possible  at  extreme 
altitudes.  This  is  not  true  because  different  conditions  are  applicable 
at  high  altitudes.  As  the  pressure  is  reduced,  the  surface  molecules 
see  fewer  and  fewer  air  molecules.  Consequently,  the  region  over  the 
surface  becomes  less  croMed  until  a  situation  is  reached  where 
essentially  all  the  molecules  which  leave  the  surface  escape.  In  this 
case  the  sublimation  rate  will  no  longer  be  limited  by  diffusion,  but 
will  tend  toward  the  absolute  rate  of  sublimation. 

w 


(27) 


Altitude  (ft) 


g.36  —  Subliming  systems:  wall  temperature  versus  mass-transfer  parameter 
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o  gravitational  constant 
a  *>  condensation  coefficient 


All  the  qmntltieB  here  may  be  calculated  vflth  the  exception  of  a, 
the  condensation  coefficient.  This  is  uemlly  determined  experimentally 
and  it  has  been  foimd  that  it  la  relatively  small  in  the  case  of  polar 
confounds  and  close  to  xinlty  for  non-polar  materials.  The  absolute 
tenperature  of  the  surface  will  also  be  an  Important  factor  for  it  will 
determine  if  a  truly  "clean”  surface  exists.  For  T  «  1000\  we  may 
assume  that  the  absorbed  layer  la  no  longer  present  at  extremely  lew 
pressures.  Therefore  the  maximiim  blowing  rate  we  may  expect  is: 


However,  in  this  case  of  rarefied  flow,  the  heating  rate  is  no 
longer  affected  by  the  efflux  of  vapor  from  the  surface.  That  is,  the 
heating  rate  la  Independent  of  the  sublimation  rate.  The  surface 
temperature  and  consequently  the  sublimation  rate,  however,  adjust 
themselves  so  as  to  balance  the  heat  input  by  heat  absorption  throu^ 
sublimation.  That  is: 


Cam  AH 
^  X  B 


(29) 


where  q  is  determined  by  flight  velocity  and  altitude  and  vehicle 
geometry. 
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Appendix  A 

ADDITIONAL  VERIFICATION  OF  THE  DIMENSIONLESS  PRESENTATIONS 

In  this  section  the  reduced  skin  frlctlen  ej,/c_  and  heat  transfer 

r'  ro 

q/q^  are  plotted  as  functions  of  the  mass -transfer  parameter  for  all  of 
the  Injected  materials#  (Figs.  37# 38*)  ^t  may  he  seen  that  the  aunnary 
curve  given  in  the  main  body  of  the  report  accurately  represents  the 
ealotilated  performance. 

Furthermore 4  the  Stanton  nimiber  and  recovery  factors  are  given  In 
detail.  (Figs.  39'*^.)  In  the  latter  case^  It  Is  apparent  that 
discrepancies  exist  In  the  case  of  bydrogen-alr  and  hellum-alr,  probably 
due  to  the  uncertainty  in  the  physical  property  values  under  such 
extreme  teiqperature  conditions. 
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g.  38  —  Effect  of  mass 
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on  heat  transfer 


Hydrogen-air 
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Fig.  40  —  Effect  of  mass 


1-2516 

85 


sr  on  Stanton  number 
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Fig.  42— Effect  of  mass  transferor!  recovery  factor 
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Appendix  B 

BEUtfTlON  BETWEEN  WALL  CONCENTRATION  AND  MASS-TRANSFER  PARAMETER 

Zn  the  treatment  of  a  eubllmlng  eurfacei.  It  is  important  to  knovr 
the  relationship  between  the  vail  mass  fraction  and  the  dimensionless 
mass-transfer  rate.  A  generalized  presentation  of  the  wall  mass 
fraction  for  each  gas  appears  possible  with  the  aid  of  the  Baron 
mriable  (Py\/Pg%)  and  these  are  shown  for  some  five 

gases  in  Figs.  l('3'-46.  A  sxsmnarizlng  curve  is  also  available  in  Fig,  46. 


Fig.43-Woll  mass  fraction  versus 
dimensionless  mass  transfer 


Woll  mass  fraction  of  injection  coolant 
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Fig.44-Wall  mass  fraction  versus 
dimensionless  mass  transfer 


Wall  mass  fraction  of  injection  coolant. 
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Fig. 45  — Wall  mass  fraction  versus 
dimensionless  mass  transfer 


Wall  mass  fraction  of  injection  coolant 
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Appendix  C 

ON  THE  USE  OF  MAS3-TOANSFER  RESULTS  FOR  PREDICTION  OF  MM  TOANSFER 

It  Is  rather  oonnnon  practice  to  obtain  experimentally  nase -transfer 
Stanton  numbers,  C^,  using  materials  which  sublime  or  evapc=>rate,  such  as 
naptheilene  or  water,  and  to  Interpret  these  results  dlreotlLyas  solid 
wall  heat-transfer  Stanton  number,  (in  some  instsmces  oa.  aoxn^eetlon 
factor  involving  the  Lewis  number  Is  applied).  However,  ecanslderahle 
caution  must  be  exercised  In  this  procedure  as  is  obvious  f  ran  a  com¬ 
parison  of  Fig.  47  with  Fig.  l8  of  the  main  text.  Figure  MS  shews  the 
heat-transfer  Stanton  numbers  in  the  form  of  where  Is  the 

heat-transfer  Stanton  number  in  the  presence  of  mass  traneSer,  while  C,, 

Bo 

is  the  Stanton  number  for  a  solid  wall  exposed  to  the  sane  free  stream 
conditions  (l.e.,  °Ho  is  the  limiting  value  of  as  the  idbibs -transfer 
rate  goes  to  zero).  Figure  4?  is  a  similar  curve  for  the cxbbb -transfer 
Stanton  number,  again  normalized  with  respect  to  the  vealue  at  the 
zero  mass-transfer  condition.  For  water  vapor-air  and  COg— air  mixtures 
the  limiting  Stanton  values,  and  Cj^,  are  in  good  agreeeiMnt  (l.e., 
®Ho  "  ®Mo^  since  does  not  depart  appreciably  freem  unity  over 

a  range  of  mass -transfer  rates,  it  would  appear  that  the  nwasured  mass- 
transfer  Stdnton  nuniber,  0^,  for  these  two  systems  could  laadeed  be 
Interpreted  as  solid  wall  heat-tremsfer  Stanton  nuniber,  Ojj^,  However, 
for  the  other  gas  mlxtvires,  it  must  be  concluded  that  cossrdderable  error 
would  arise  in  the  use  of  this  procedure. 


Dimensionless  moss  transfer, 


/’w 


p. 


Fig. 47—  Effect  of  moss-transfer  rote  on  the 
moss-transfer  Stanton  number 
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